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Abstract 

The Casimir force is calculated between plates with thin metallic coating. Thin films are de- 
scribed with spatially dispersive (nonlocal) dielectric functions. For thin films the nonlocal effects 
are more relevant than for half-spaces. However, it is shown that even for film thickness smaller 
than the mean free path for electrons, the difference between local and nonlocal calculations of 
the Casimir force is of the order of a few tenths of a percent. Thus the local description of thin 
metallic films is adequate within the current experimental precision and range of separations. 



1raul@fisica.unam.mxl 

i V.B.Svetovoy@el.utwente.nl On leave from Yaroslavl University, Yaroslavl, Russia 



1 



I. INTRODUCTION 

The Casimir force between uncharged metallic plates [l| (see also reviews DDDDDD 
attracted considerable attention in the past years. The force was measured in a number of 
experiments with a hi g h p recision using different techniques and geometric configurations 
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On the other hand, the potential applications of the force in 



micro and nanomechanics is still largely unexplored. Actuation and nonlinear behavior of a 



and the importance of 
. Due to technological 



mechanical oscillator with the Casimir force were demonstrated 
the force in adhesion and stiction has also been discussed jl?! [l8|, [it 
reasons thin coating layers or multilayered structures are often in use in micromechanical 
devices. The main question to be addressed in this paper is how important are the nonlocal 
effects when the film thickness is smaller than the mean free path of the electrons. 

For the first time the problem of a thin metallic layer on top of another metal appeared 
in connection with the first atomic force microscope (AFM) experiments P Jlo| . In these 
experiments a relatively thick Al layer was covered with Au/Pd film of 20 nm jj.9] or 8 nm 
thick to prevent aluminum oxidation. Because the film was thin enough to be transparent 
for the light with a characteristic frequency u c h = c/2a, where a is the distance between 
the bodies, it was concluded that Au/Pd layer did not influence on the force 20|. In actual 
calculations [2l( the thin Au/Pd film was changed by vacuum. This approach was criticized 



22j on the basis that according to the Lifshitz formula (23, UM ^ ne force depends on the 



dielectric function e at imaginary frequencies u = i(. Kramers-Kronig relation shows 
that at ( c h = c/2a low real frequencies uj <C c/2a give significant contribution to £(i( c h)- 
At low frequencies Au/Pd film is not transparent and it should be taken into account. It 
was demonstrated that, indeed, even 8 nm thick film gave significant contribution to the 



force. The calculation in Ref. 



251 ] supported this conclusion but the authors speculated 



that nonlocal effects due to small thickness of the film (smaller than the mean free path for 
electrons) allowed one to consider the film as transparent. 

The question arose again in connection with a recent experiment j^J, where the force 
was measured between a plate and sphere covered with 10 nm or 200 nm Pd film. For thin 
film the expected reduction of the force was clearly observed. It was indicated that spatial 
dispersion might be im por tant for calculation of the force in the case of thin film. Also, 
Bostrom and Sernelius [2J| pointed out to the need of detailed studies of nonlocal effects, 
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while studying the retarded van der Waals force between thin metallic films within a local 
approximation. 

There have been several works dealing with the problem of nonlocality in the Casimir 
force between half spaces. Katz JscJ was the first to point out the need of a quantitative 
study and in his work only a rough estimate of how spatial dispersion affected dispersive 
forces was given. Heindricks jjyj was able to derive Lifshitz formula in an approximate way 
to include nonlocal effects. Similarly, Dubrava using a phenomenological approach described 
the Casimir attraction between thin films 3^|. More recently, based on the formalism of 
nonlocal optics, the effects for thick metallic layers have been considered. Propagation of 



bulk plasmons 



3, 24, [35 



35l ] and electromagnetic response in the region of anomalous disper- 



sion [36| were taken into account, showing that the spatial dispersion does not contribute 



significantly to the Casimir force. The method developed in Ref. 



is very general and 



can be used for the analysis of all nonlocal effects including those arising in thin films. 



II. FORMALISM 



The Casimir force between two plates separated by a vacuum gap a at a temperature T 
is given by the Lifshitz formula j3, E | . The force is expressed via the reflection coefficients 
Ri and R 2 of the plate 1 and 2, respectively, in the following way: 



Fpp (a) 



k B T 
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' / dqqk {(R^R^ 1 exp (2ak ) - l) 1 + (R^R^ exp (2ak Q ) - l) 



n=0 



-1 



(1) 



where subscripts s and p denote the polarization states, q is the wave vector along the plates, 
q = |q|, and k is the normal component of the wave vector defined as 



q 2 . 



ko=VQ/c^T 

In Eq. (HJ) the sum is calculated over the Matsubara frequencies 



(2) 



. 2TTk B T 



(3) 



The reflection coefficients R\ and R2 are different for s and p polarizations and are functions 
of q and imaginary frequencies ( n . They comprise material properties of the plates and for 



this reason we start our analysis from the reflection coefficients. 



A. Local case 

To set our notation, we first study briefly the known local case, when the optical response 
depends only on frequency. We start the analysis from a thin film of thickness h on a 
substrate. It will be assumed here that the film is continuous. 

For a film on an infinitely thick substrate the problem is rather simple. The Maxwell 
equations are solved with the boundary conditions which are the continuity of the tangential 
components of electric and magnetic fields on both boundaries of the film. The problem can 
be solved at real frequencies and then analytically continued to the imaginary axis. In the 
local limit the film and substrate are described by their local dielectric functions which will 
be denoted as E\ (u) and e 2 (^)- in general, the indexes marking the layers will increase from 
top to bottom of the plate. The dielectric function of vacuum will be taken as Eq (uj) = 1. 
The reflection coefficients in our case are well known in optics. At imaginary frequencies 
they are [37| : 

R _ 7"oi - ?"2iexp (-gfcife) 
1 - ^01^21 exp (—2k\h) ' 

where r m i are the reflection coefficients from the boundary between media I and m. These 
coefficients depend on the polarization, s or p, and are defined as 

m "'I ^p E-lkm £"nv*>l i^i) 

k m 4~ k\ £ik m + s m ki 

where k m is the normal component of the wave vector in the medium m: 



C 2 

k m = d e m (i() — + q 2 . (6) 

It is easy to check that for h — > 00 (thick film) R — > roi and in the opposite limit h — > the 
reflection coefficient coincides with that for the substrate: R —> r 02 . 

To understand the variation of the Casimir force with the film thickness, we first study 
the behavior of the reflection coefficients. For a qualitative analysis it will be assumed that 
a metal, film or substrate, can be described with the Drude dielectric function 
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where cj p and u T are the Drude parameters which are different for each layer. For thin 
films, u T is a function of the film thickness. This dependence appears because, in addition 
to the internal scattering processes, for thin films scattering from the surfaces is important. 
These processes are independent of each other and the relaxation time in the Drude model is 
uo T = uj^ ulk + u^ ur f(h). This effect becomes important when the thickness is smaller than the 
mean f ree p ath for electron. Dependence on h of u s T ur ^ is explained by the Fuchs-Sondheimer 
theory [ssi Isoj]. When h is much smaller than the mean free path, this dependence is given 
by 

< urf (h) = l(l- P )^, (8) 

where vp is the Fermi velocity and an electron has probability p of being specularly reflected 
from the surface. As one can see from Eq. (jHJ) only diffusely reflected electrons contribute to 
uj s T ur f(h). Experimental results concerning the specularity are far from unique. Very different 
values of p in the range < p < 1 were used to explain the experimental results 40] . In this 
paper we investigate the nonlocal effects for specular reflection of electrons on the surface 
and do not include in the consideration /i-dependence of the relaxation frequency. But in 
any case our results are not very sensitive to the exact value of tu T . 

Consider first the system consisting of S1O2 substrate with 62 = 4 and Au film on top 
of it with the parameters u p = 9.0 eV, uj t = 0.035 eV [41]. It is convenient to introduce 
dimensionless variables and parameters as follows: 

n C n c 1 uj p h 

^ = — , Q = — , 7 = — , H = . (9) 

UJp UJp LOp C 

The reflection coefficient for s-polarization as a function of the dimensionless frequency Q 
is shown in Fig. ^ The dashed curve corresponds to semi-infinite metal h — ► 00. It was 
calculated with Q = 0.1. This value is taken for the characteristic wave number q ~ I /2a 
at a ~ 100 nm. The solid lines marked as 1, 2, and 3 correspond to the dimensionless 
thickness H = 0.3, 1, and 3, respectively. Note that H = 1 gives the film thickness h equal 
to the penetration depth 5 = c/uj p w 22 nm (Au). One can see that R s decreases fast with 
the thickness. When Q increases the film also becomes more transparent for s-polarization. 
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FIG. 1: Reflection coefficient for s-polarization as a function of SI in the local case (metal film on 
the dielectric substrate). All the results are presented for Q = 0.1. Infinitely thick film is given by 
the dashed line. The curves 1, 2, and 3 correspond to the dimensionless thickness H = 0.3, 1, and 
3, respectively. 

The other distinctive feature is that R s is going to zero in the limit SI — > 0. In this limit 
s-polarized field degenerates to pure magnetic field, which penetrate freely via the metallic 
film. 

The reflection coefficient for p-polarization shows a different behavior as one can see in 
Fig. |21 The dashed line represents the thick film and the solid lines marked as 1 and 2 
correspond to H — 1 and 0.1, respectively. Variation of R p with the film thickness is not 
very significant. The reason for this is the effective screening of the E z component even by 
a very thin metallic layer. An important conclusion can be drawn from this simple fact. 
The film thickness affects mostly the contribution of s-polarization, but the part of the force 
connected with p-polarization is changed weakly in the local case. 

Consider now the effect of a thin film on top of a thick metallic layer. It will be assumed 
that both metals can be described by the Drude dielectric functions e\ and e<i which 
differ from each other only by the values of parameters u>i P and u>i T (i = 1, 2). It is clear that 
in dependence on the film thickness the reflection coefficients will be in between the lines 
describing metal 1 (h — > oo) or metal 2 (h — > 0). In Fig. El we present the case when the 
top layer is better reflector than the bottom one. The dotted line gives r 02 and the dashed 
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FIG. 2: Reflection coefficient for p-polarization as a function of Q in the local case (metal film on 
the dielectric substrate). Infinitely thick layer is shown by the dashed line. The solid lines marked 
as 1 and 2 correspond to H = 1 and 0.1, respectively. 

line represents tqi- The results for the film with thickness H = 1 and 0.1 are marked as 1 
and 2, respectively. In our calculations, the ratios uji p /uj2 P = 2 and oj\ t /uj2t = 1 where used 
and dimensionless parameters (jHJ) were denned relative to the parameters of the top layer 1. 
The relaxation frequencies, a>j r , influence mostly on low frequency behavior of R s . They are 
not very important for the Casimir force because the main contribution in the force comes 
from the imaginary frequencies Q ~ c/2au p ^> 7 where u T does not play significant role. 
The reflection coefficient, R p , for p-polarization is shown in Fig. 0] The curves 1 and 2 
correspond H = 1 and H = 0.1, respectively. Again one can conclude that the top layer is 
more important for s than for p-polarization. 



B. Nonlocal case 



For propagating photons the reflectivity of thin films in the nonlocal case has been ana- 
lyzed in Ref. |42j. It was assumed that electrons are reflected specularly on both boundaries 
of the film. Let us consider first s-polarization. Similar to the case of a semi-infinite metal 
(43 1 the tangential component of the electric field is considered as even on each boundary: 
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FIG. 3: Reflection coefficient for s-polarization as a function of SI in the local case (metal film on 
the metallic substrate). The dotted line represents the substrate, H = 0, the dashed line represents 
the thick top layer, H — > oo. The curves 1 and 2 correspond to H = 1 and H = 0.1, respectively. 

1 

0.8 
0.6 

Q. 

0.4 
0.2 

10~ 5 10~ 3 10"' 10 1 

Q. 

FIG. 4: The same as Fig. Olbut for p-polarization 

E y (rah - z) = +E y {rah + z) , (10) 

where z is the direction normal to the film surface, ra is an arbitrary integer, and the plane 
of incidence was chosen to be x — z. The Maxwell equations and Eq. (flT)|) demand for the 
magnetic field on the boundaries the following conditions: 

8 




H x (mh — z) = —H x (mh + z) , H z (mh — z) = +H Z (mh + z) . (11) 

Formally the conditions (JTUJ), (jH| continue the film of finite thickness to the infinite layer. 
These conditions mean that the fields can be considered as periodic with period 2h, and 
they can be expanded in a Fourier series. 

In the nonlocal case the material is characterized by the impedance instead of local 
dielectric function. The impedance is defined as the ratio of tangential components of 
electric and magnetic fields just below the surface. For s and p-polarizations the impedances 
of metallic film were found in Ref. 4^ with the method which is direct generalization of 



the method used for semi-infinite layer 
one can define on each of them: 



431 ] . The film has two surfaces and the impedances 



Ey 



z=5,h—S 



Ex 



(12) 



z=5,h— S 



where S — > 0. It was noted 42] that instead of impedances ()T2j) one can use a different couple 
for each polarization which can be easy calculated. These new impedances were introduced 
as the ratio of the fields even or odd relative to the film center z = h/2. Even or odd fields 
will be marked by the superscripts (1) or (2), respectively. The new impedances 



Z(1.2) 



E„ 



(1,2) 



H 



(1,2) 



Z (l,2) 
P 



E.. 



(1.2) 



(1,2) 



are the same on both boundaries of the film because of the symmetry conditions 



(13) 



(5) = 43 (h-S), Efl (5) = -Eg} (h - 5) 

and similarly for the magnetic field. 

Explicit expressions for these impedances were found in Ref. 42j: 



(14) 



Ck uneven) k ) ~ (f ) " ^ 



z{ 2) = .2U y, 1 

p ch ^ k 2 

n=(odd, even) 



I h ) 



(15) 



(16) 



where for even, (1), or odd, (2), fields the sum has to be calculated over n = 2m + 1 
or n = 2m, respectively. The transverse dielectric function e t (uj, k) contributes to Z s . It 



describes the response of the material on the electric field transverse to the wave vector k. 
In case of the p-polarization z-component of electric field creates a nonzero charge density 
in the metal producing the longitudinal field inside of metal. That is why Z p depends also 
on the longitudinal dielectric function Ei (u,k). In general, these functions are nonlocal, so 
they depend on both uj and k. The absolute value of the wave vector k in Eqs. (|15|) . is 



Let us consider now the reflection and transmission coefficients of the film on a substrate. 
Note that in Ref. only a free standing film was considered. To find these coefficients 
one has to match the tangential components of the electric and magnetic fields outside and 
inside of the film. We assume for simplicity that the substrate can be described by a local 
dielectric function or equivalently by local impedances. This assumption is justified by the 
investigation of nonlocal effects at imaginary frequencies for semi-infinite metals [36]. It was 
demonstrated that in contrast with the real frequencies the nonlocal effect (anomalous skin 
effect) brings only minor influence on the reflection coefficients. Matching the electric field 
on both sides of the film for s-polarization one gets 

Ef (1 + R a ) = (6) + 4 2) (6) , 
E° y t s e ik * h = EP (5) - Ei 2) (5) , 

where is the incident field, t s is the transmission coefficient and the symmetry conditions 
(114}) were taken into account. Similar equations are true for the magnetic field: 

H° x (1-R s )=hP(5)+hP(8), 

Eqs. and (JT5j) can be solved for R s and t s using the impedance definition (JT3J). As the 
result the reflection coefficient can be presented in the form: 



(jff - (g + 2,2) + (g - z, ) (z<j> + z*\ 

R ' " (Z,<;> + zj) (ZS> + zj) + (Zff + zj) + Z. X ' <20) 



Here we introduced the following notations : Z sl ' are the nonlocal impedances of the film 
given by Eq. (jT3j) . Z s2 is the local impedance of the substrate defined as 
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and 



Z s0 = (22) 
ck 

is the "impedance" of the plane wave defined as the ratio of electric and magnetic fields in 
the wave. The formula ()20j) for R s cannot be presented in the same form (0J) as in the local 
case. This is because we used the impedances ()13|) instead of that given by Eq. (fT2~|). As we 
will see both Eqs. (|3J) and (J2U|) coincide in the local limit. 

In the same way one can find the reflection coefficient for p -polarization, R p . In this case 
the equations similar to (|TH|) . (JTHJ) with the interchange x y will be true, the impedance 
of the plane wave is defined as 

Z p0 = — = ^— , (23) 
and the local impedance of the substrate is 

ck 

^p2 — — 

id 

The final expression for R„ is 



uje 2 (oj) 



(24) 



{z$ + z p0 ) (z® + z p2 ) + (z^ + z p0 ) (z« + z p2 ) ' 

It differs from Eq. (J2~U|) only by the general sign and the change s — > p. 

If the substrate is changed by vacuum, Z a2 — >■ Z^o (« = s, p), we reproduce the reflection 



coefficient found in Ref. 



4J: 



R a = \ (ri 1} + rf) , a = s, p, (26) 

where the "partial" reflection coefficients are connected with the impedances by the usual 
relations 



7 7(1,2) 7 7 (i,2) 

(1,2) ^0 ~ z * r (i,2) ^po / 97 \ 
s " 7 ,7(1,2)' 'p 7 , 7 (i,2)- 
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In the local limit both the transverse St and longitudinal Si dielectric functions coincide with 
the local function: e t (o>, k) — > E\ (u, k) — > e% (u). In this case the sums in Eqs. (JT3J), (JTfiJl 
can be found explicitly. For example, for s-polarization one has 

ls cki 2 ' ls cki 2 v ; 

Substituting it in Eq. ()20j) one can check that the reflection coefficient for the local case 
given by Eq. (J3J) is reproduced. 

All the equations above were written for real frequencies. Transition to imaginary frequen- 
cies, which are the main point of our interest, can be done by a simple analytic continuation. 
To get the nonlocal effects in the reflection coefficients, we have to fix the nonlocal dielectric 
functions. At imaginary frequencies in the Boltzmann approximation they are given by the 
relations [36[ 

,^ s fi ( v ) n , s 3 v — arctan-y . 

£ ^ = 1 + nWWY / ^ ) = ^+^-arctan,r W 

£t v) = 1 + n{n + 7 y ft {v) = 2^ + ^ + v ^ arctany ] ' ( 3 °) 



~c Q + 7 ' (31) 
where is the Fermi velocity. The dimensionless variables Q have been introduced in 
Eqs. (|29|) -([3"T |) . In addition, we have neglected in Eqs. (|2*5j) -(|3Uj) the contribution due to 
the interband transitions. 

The reflection coefficients in the nonlocal case were calculated numerically. In Fig. Elthe 
difference between local and nonlocal coefficients R s is shown for Au film on top of S1O2 
substrate. The dashed curve corresponds to very thick film, H — ► 00. The solid lines marked 
as 1 and 2 are presented for H = 1 and 0.1, respectively. As before, H = 1 corresponds to the 
penetration depth of Au (5 = 22 nm). The thick film clearly demonstrates the anomalous 
skin effect at Q ~ 7, although the magnitude of the effect is small as was already noted in 
Ref. |36j]. Even this small effect decreases with the film thickness as the curves 1 and 2 show. 
The nonlocal effect increases with Q but it is smaller than 1% even for Q — 1. It should 
be noted that the Boltzmann approximation is good while Q < 1, but when Q approaching 
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FIG. 5: Difference between the reflection coefficients for s-polarization in the local and nonlocal 
cases (metallic film on the dielectric substrate). The dashed curve was calculated for very thick 
film. The solid curves are presented for H = 1 (1) and H = 0.1 (2). 

1 the reflection coefficient itself becomes small and there is no sense to keep the nonlocal 
correction in this range. Similar result was found for the film on top of a metallic substrate. 
One can conclude that for s-polarization the nonlocal effect in the reflection coefficient is 
very small and can be neglected in calculation of the Casimir force. 

The situation for p-polarization is shown in Fig. |U]for the film on top of metallic substrate. 
As in the local case the substrate was chosen to have the plasma frequency 2 times smaller 
than that for the film. The lower curve corresponds to Q = 0.1 and H = 1. The upper series 
of curves is given for Q = 0.5. As one can see, the nonlocal effect manifests itself in a wider 
frequency range and does not disappear even for zero frequency. The latter is the result 
of Thomas- Fermi screening as was explained in Ref. [^. The effect is still small but the 
nonlocal contribution in the Casimir force will be larger than that for s-polarization. This 
is because the nonlocal effect is the largest at frequencies which give the main contribution 
in the Casimir force. 
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FIG. 6: Difference between the reflection coefficients for p-polarization in the local and nonlocal 
cases (metal film on the metallic substrate). The lower line is given for Q = 0.1, H = 1. The upper 
series is for Q = 0.5. The dashed curve presents very thick film. The solid curves are given for 
H = 1 (1) and H = 0.1 (2). 

III. EFFECTS OF SPATIAL DISPERSION ON THE CASIMIR FORCE 

To quantify the effect of spatial dispersion on the Casimir force, we calculate the percent 
difference between the local case and nonlocal case (A% = \(Fi OC(d — F non i oca i) / Fi ocal \) , as a 
function of separation. 

First we consider the case of free standing metallic films. The system is similar to that 
considered by Bostrom and Sernelius j^J. The percent difference A% as a function of 
separation is presented in Fig. Q for three different thicknesses. The results for the thick 
ilm h = 100 nm coincides with the results obtained for half spaces in our previous work 
3(| . As the thickness decreases the nonlocal effects become more relevant. Thin films 
have a more complicated nonlocal response than half spaces. For p-polarized waves, surface 
plasmons on each side of the film can interfere [? ], creating standing waves that will 
increase the electromagnetic absorption of the field that will decrease the Casimir force. 
These resonance conditions are evident from Eq. (fTTj) where k z = nn/L. 

The force is not affected significantly when the thin films are on substrates. In Figure[B]we 
have plotted the percent difference between two thin Au films, each deposited on a dielectric 
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FIG. 7: Percent difference A% between the local and nonlocal Casimir force between Au free 
standing films. The values for the thick film coincide with those obtained for half-spaces 
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Separation(nm) 



FIG. 8: Percent difference A% between the local and nonlocal Casimir force between two Au films 
deposited on a dielectric substrate. 

substrate. Again, we assumed e = 4 for the dielectric, just as an illustrative example of the 
effect of substrate. The substrates reduce slightly the value of A% for both curves shown, 
with the obvious limit that when the substrate has the same dielectric function as the film, 
we recover the results for the force between half-spaces. This means that the effect of the 
substrate is to allow energy transfer out of the thin film into the substrate. 

The difference between the local and nonlocal cases can be reduced in a system consisting 
of Au half space and Au coated substrate. Again, we took a dielectric (e = 4). The effect 
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FIG. 9: Percent difference A% between the local and nonlocal Casimir force between Au half space 
and Au coated dielectric. 

TABLE I: The magnitude of the nonlocal effect for different substrates at a fixed separation of 50 
nm and a film thickness of 2 nm. 



substrate 




A% 


dielectric, e = 


= 4 


0.34 


metal, = 


2 


0.37 


metal, = 


0.5 


0.44 


no substrate, e 


= 1 


0.44 



of spatial dispersion reduces significantly as compared to the cases treated in Figs. [7| and |S1 
This shows that the most important part of the spatial dispersion effect come from the thin 
films. If in current experiments the separation can go down to 50 nm, in the system shown 
in Fig. EH the nonlocal correction is of the order of 0.34%. 

The result holds for different substrates. This is shown in Table 1, where we presented the 
percent difference between the local and nonlocal forces for Au film deposited on different 
substrates. All data are given for a separation of 50 nm. As before, uj p i/uj P 2 is the ratio of the 
Au plasma frequency to that of the metallic substrate, assuming the damping factor remains 
the same. The case e = 1 corresponds to the free standing Au thin film (no substrate). 
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IV. CONCLUSIONS 



The role of thin metallic coatings in the calculation of Casimir forces has been studied 
taking into account spatial dispersion. The description of the nonlocal response of thin films 
is based on the Kliewer and Fuchs formalism that imposes a symmetrical behavior of the 
fields inside thin films. The study of the reflectivities shows that the main contribution 
to the nonlocal effect comes from p-polarized light that excites normal modes within the 
material. At very small separations, the effects can be appreciable but at best a percent 
difference of 7% is found. However, for typical experimental setups and separations the 
percent difference between the local and nonlocal case is of the order of 0.4%, that can 
be regarded as negligible within current experimental precisions and the local description 
is good enough. The effect of thin films within a local approximation has been measured 
recently by Lisanti et al. 



34, 



35 



wc 



Along with the previous works on nonlocal effects between half-spaces |33 
can generally conclude that these effects will be difficult to detect at the current experimental 
precision. Our results indicate a decrease in the force due to spatial dispersion. However for 
half-spaces within a jellium model it has been shown [35j that the force can increase due to 
nonlocal effects because of decrease in the separation of the optical surfaces that might not 
coincide with the physical surface. 
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